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ABSTRACT 

We reexamine the general synchrotron model for GRBs' prompt emission and de- 
termine the regime in the parameter phase space in which it is viable. We characterize 
a typical GRB pulse in terms of its peak energy, peak flux and duration and use the 
latest Fermi observations to constrain the high energy part of the spectrum. We solve 
for the intrinsic parameters at the emission region and find the possible parameter 
phase space for synchrotron emission. Our approach is general and it does not de- 
pend on a specific energy dissipation mechanism. Reasonable synchrotron solutions 
are found with energy ratios of 10~^ < ts/te < 10, bulk Lorentz factor values of 
300 < r < 3000, typical electrons' Lorentz factor values of 3 x 10^ < 7e < 10^ and 
emission radii of the order lO^^cm < R < lO^^cm. Most remarkable among those 
are the rather large values of the emission radius and the electron's Lorentz factor. We 
find that soft (with peak energy less than lOOKeV) but luminous (isotropic luminosity 
of 1.5 X 10^^) pulses are inefficient. This may explain the lack of strong soft bursts. 
In cases when most of the energy is carried out by the kinetic energy of the flow, such 
as in the internal shocks, the synchrotron solution requires that only a small fraction 
of the electrons are accelerated to relativistic velocities by the shocks. We show that 
future observations of very high energy photons from GRBs by CTA, could possibly 
determine all parameters of the synchrotron model or rule it out altogether. 

Subject headings: gamma rays: bursts, theory, method: analytical - radiation mecha- 
nisms: non thermal 



1. Introduction 



GRBs' non thermal spectrum has lead to the suggestion that both the prompt and the afterglow 
are produced by synchrotron emission of relativistic electrons ( |Katz||199"4| |Rees and Meszaros 
1994[ |Sari et al.]|1996[ |1998| ). There is reasonable agreement between the predictions ( |Sari et alT 
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1998p of the synchrotron model and afterglow observations ([Wijers and Galamal|1999[ [Granot et 



al.|1999[[2002t|Panaitescu and Kumar(200T||Nousek et al.|2006t [Zhang et al.|2006l ). However, the 
situation concerning the prompt emission is more complicated. 

The observations of many bursts whose lower spectral slope is steeper than the "synchrotron 
hne of death" ( [Crider et"aL]|1997[ [Preece et"aL]|1998[ |2002| ) pose a serious problem to synchrotron 



emission. This has motivated considerations of photospheric thermal emission (Eichler & Levinson 
2000t|Meszaros fc Rees|2000t|Peer et al.|2006t[Rees fc Meszarosl2005||Giannios|2006i [Thompson 
et al. 2007[ Ryde and Peer |2009| ), in which thermal energy stored in the bulk is radiated in the 
prompt phase at the Thomson photosphere and the high energy tail is produced by inverse Comp- 
ton. While this model yields a consistent low energy slope it has its own share of problems. Most 
notably Vurm et aL] (2012[) have recently shown that unless the outflow is moving very slowly no 
known mechanism can produce the needed soft photons. Furthermore, it is not clear how the high 
energy GeV emission can be produced within an optically thick regime ( Vurm et"aL][2012a ). On 
the other hand, observations of some bursts, e.g. GRB 0809 16C in which there is a strong upper 
limit on the thermal component ( [Zhang et alT][2009l ), or GRBs: 100724B, 110721A, 120323A, 
where a thermal component was possibly detected but with only a small fraction (5 — 10%) of 



the total energy in the thermal component ( jGuiriec et al. 2011[ 2012, Axels son et al.|[20T2 ) and 
other bursts whose spectral slopes is consistent with synchrotron suggest that for some bursts syn- 
chrotron might be a viable model. In fact, in all cases that involve Poynting flux outflows that are 
Poynting flux dominated at the emitting region synchrotron is inevitable and possibly dominant 
( |Beniamini k. Piran][2013[ ). As such it is worthwhile to explore the conditions needed to generate 
the observed prompt emission via synchrotrorj^ . 

To examine the conditions for synchrotron to operate, we consider a general model that fol- 
lows the spirit of Kumar & McMahon ( 2008[ ). We don't make any specific assumptions concerning 
the details of the energy dissipation process or the particle acceleration mechanism. Instead we as- 
sume that relativistic electrons and magnetic fields are at place within the emitting region, that 
is moving relativistically towards the observer. We examine the parameter space in which these 
electrons emit the observed radiation via synchrotron while satisfying all known limits on the ob- 
served prompt emission. We use a simple "single-zone" calculations that do not take into account 
the blending of radiation from different angles or from different emitting regions (reflecting matter 
moving at different Lorentz factors), but for the purpose of obtaining the rough range of relevant 
parameters this is sufficient. 



'Note that the synchrotron-self Compton (Waxman 



1997 



high energy LAT observations (Piran et al. 2009 Zou et a^ 



Ghisellini and Celotti 



1999 I have been ruled out by the 



2009 | l, while external inverse Compton ( |Shemi|[T994 



|Brainerd[[1994[ |Shaviv and Dar|[T995[ |Lazzati et al.[[2003] l require a strong enough external source of soft photons, 
which is not available in GRBs. 
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We characterize the conditions within the emitting region by six parameters: the co-moving 
magnetic field strength, B, the number of relativistic emitters, Ng, the ratio between the magnetic 
energy and the internal energy of the electrons, e = e^/ee, the bulk Lorentz factor of the emitting 
region, F, the minimal electrons' Lorentz factor in the source frame, 7^ and the ratio between 
the shell crossing time and the angular timescale, k. We then characterize a single, "typical", 
GRB pulse, which is the building block of the GRB lightcurve, by three basic quantities: the peak 
(sub-MeV) frequency, the peak flux and the duration. We compare the first two, the observed 
peak (sub-MeV) frequency and the peak flux with the predictions of the synchrotron model and 
we use the observed duration of the pulse to limit the angular time scale ( |Sari and Piran 1997a| ). 
The synchrotron solution is incomplete without a determination of the accompanying synchrotron 
self-Compton (SSC) emission, which influences the efficiency of the synchrotron emission. Fur- 
thermore, the observations of the high energy SSC emission poses additional constraints on the 
emitting region. We describe, therefore, a self-consistent synchrotron self-Compton solution. 



Three additional constraints should be taken into account (see also Daigne et al.|201 First, 
energy considerations pose a strong lower limit on the efficiency. The energy of the observed sub- 
MeV flux, is huge and already highly constraining astrophysical models. Thus, to avoid an "energy 
crisis", the emission process must be efficient and emit a significant fraction of the available energy 
in the sub-MeV band. Second, the emission region cannot be optically thick. Additionally, the 
LAT band (100 MeV- 300 GeV) GeV component is significantly weaker in most GRBs than the 
sub-MeV peak ( [Granot et aL][2009| |20T0t [Beniamini eTaL]|20TT| [Guetta et air][20TT| [Ackermann 
et al.||2012l ). We combine these three conditions and constrain the possible phase space in which 
synchrotron can produce the prompt emission. 

The paper is organized as follows. We describe in ^the basic concepts of the model, dis- 



cussing the observations in ^: 2. 1 and the structure of the emitting region and the parameter phase 



space in ^:2.2 In Owe describe the synchrotron equations (^3.1 1, their solutions (^3.2), the con 



straints from the accompanying SSC (^3.3), the implications of the results {[3.4) and energy and 
efficiency constraints ( §3.5[ ). We continue with a discussion of specific issues concerning the (i) 
expected spectral slopes and possibilities to alleviate the issue of the "synchrotron line of death" 



([3.6) and (ii) the narrow Epeak distribution (^3.7 1. In ^|4| we deviate from the general philosophy 
of the paper and we discuss synchrotron emission within the context of the popular internal shocks 
framework. We turn, in ^|5| to the implications of the pair creation limits from observed GeV emis- 
sion for the synchrotron model. In ^ we briefly discuss the possibilities of detecting GRBs with 
CTA, and their implications on the model at hand. 
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2. The Global Picture 
2.1. The observations 



Limits on the optical (Roming et al. 2006, 


Yost et al. 2007, Klotz 


2009), x-ray (O'Brien et 


al. 2006), GeV (Ando et al.|2008 Guetta et al. 


2011 


Beniamini et al. 


2011, Ackermann et al. 


2012P and TeV (Atkins et al.|2005 ; Albert et al.||2006 , . 


A.haronian et al.||2009) emission during the 



prompt phase of GRBs, demonstrate that the observed sub-MeV peak carries most of the GRB's 
energy. The only (unlikely) possibilities that have not been ruled out are either an extremely strong 
and sharp peak between 10-lOOeV or a peak at extremely high energies above the TeV range. 
Therefore, we associated the sub-MeV peak with the peak of the synchrotron emission. 

We consider the typical observations after redshift corrections, i.e. in the host galaxy frame 
which we denote as "the source frame" (as opposed to "the observer frame"). The basic ob- 
servables are: the peak frequency: hvpeak = hvpeak,obs x (1 + Z) (300KeV), the duration of 



a typical pulse, tp = tp,o6s/(l + Z) (0.5 sec) and the peak flux, F^^^^^ 



^peak ,obs 



/(I + z) 



(1.5 X 10~^^erg/cm^sec Hz). With these definitions we remove the (1 + Z) dependence from 
the equations, and the only remaining dependence on redshift is through the luminosity distance, 
relating the flux and luminosity. All three quantities vary from pulse to pulse and the values in 
brackets denote the canonical values that we use here. These three observations together with 
a typical luminosity distance, (11 = 2 x 10^® cm (or 2; = 1) correspond to an isotropic equiv- 
alent luminosity of Li^o ~ 2 x lO^^erg/sec. It has been suggested that Ep^ak and L satisfy: 
Lp = 10^^-^^ {Epeak{l + ^)/355KeV)^-^erg/sec ( [Yonetoku etaL]|2010[ ), we therefore choose Up^ak 



and Lp such that they are compatible with this relation. In [3.1 we explore the dependence of 
the synchrotron model on this choice and look at 4 representative "GRB types". An additional 
observation that we use is the limit on the total flux observed in the LAT band (30 MeV-300 GeV), 
which is at most 0.13 of the GBM (8KeV-40MeV) flux ( |Beniamini eraL][20TT| |Ando et al.|[2008 
Guetta et al.|201 m Ackermann et al.|2012D . 



2.2. The emitting region 



The general model requires a relativistic jet, traveling with a Lorentz factor T with respect to 
the host galaxy with a half opening angle of 9, at a distance R from the origin of the explosion and 
with a kinetic energy Etot in the source frame. The source can be treated as spherically symmetric 



as long as r ^ < 6 which is expected to hold during the prompt phase, in which F ^ 100 (Fenimore 
eraL][T993t [Woods and Loebl[T995t [Piran|[T995[ |Sari and Piran||T999| ). We therefore do not solve 
for 9 and use isotropic equivalent quantities throughout the paper. An upper limit on the emitting 
radius is given by the variability time scale, tp ^ tang = -R/(2cF^). We define a dimensionless 
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parameter k such that kR/G'^ is the shell's width: 



kR 



kt 



(1) 



t, 



cross 



2cr2 



ang- 



With this definition, the pulse duration is: 



tp 



{k + l)R 

2cr2 



(2) 



Clearly, for = 1 the pulse width is determined by angular spreading. 

We consider a single-zone model in which the different properties of the emitting region are 
constant throughout the zone. Notice that the single-zone approach cannot test spectral evolu- 
tion, and becomes less accurate as k increases. Initially, the emitting region consists of either an 
electron-proton or an electron-positron plasma. The numbers of protons and positrons are denoted 
by Np and Ne+ respectively. Due to electrical neutrality, the number of electrons is Np + Ne+. 
After the gamma-rays are produced, annihilating high energy photons may result in extra electrons 
and positrons. We denote the number of pairs created in this way by N pairs (typically these par- 
ticles are less energetic than the original population of emitting particles). Altogether, there are 
Ntot = Np + 2Ne+ + 2Npairs clcctrons and positrons in the flow, the number of electrons is always 
larger than the number of positrons except for the case of an initial pair dominated outflow, where 
they are equal. Therefore, to simplify the presentation, we use hereafter "electrons" instead of the 
more accurate but cumbersome "electrons and positrons". In the fast cooling regime, the typical 
synchrotron cooling time of a relativistic electron is shorter than the pulse duration and thus the 
instantaneously emitting particles are much less than the total number of emitting particles during 
one pulse (see Fig[T]). The emitting electrons are characterized by X^s < 1 which is defined to be 
the ratio between the number of instantaneously emitting electrons and the overall number of rela- 
tivistic electrons emitting during the entire pulse (A concise definition of all these particle numbers 
is given in table [T]). 



Notation 


Description 


Ntot 


The total number of electrons (and positrons) 


Npairs 


The number of pairs created by annihilating high energy photons ' 


Np 


The number of protons in the initial flow 


Ne+ 


The number of positrons in the initial flow 


Ne 


The number of relativistic electrons (and positrons) radiating during one pulse 


Xins 


Ratio between instantaneously emitting electrons (and positrons) and Nf> 



Table 1 : Different particle numbers used in the text. 

^Excluding the original pairs that may reside in the initial flow and are accounted for by A'^e^- 
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Instantaneously 
emitti ng elect rons 



uniform magnetic field — B 




kR 

thicknes — — 77 

number of relativistic electrons = 
number of electrons radiating during tp = 







Fig. 1 . — A schematic representation of the single zone model. The model assumes that the physical condi- 
tions are uniform throughout the emitting region. Due to the short cooling time by radiation (compared with 
the pulse duration) the instantaneously emitting particles reside in a thinner slab than the overall population. 



The relativistic electrons have a power-law distribution of energies: 



(3) 



(where C is a normalization constant) which holds for 7m < 7 < iMax- Based on theory, we expect 
p to be of the order p ~ 2.5 ( [Achterberg et al.|[2001[ |Bednarz et al.||1998| [Gallant fc Achterberg 
1999| ). Indeed this result has been confirmed observationally both from the GRB prompt and 
afterglow phases ( |Sari and Piran|[T997b| [Panaitescu and Kumar||2001 [). 7Max arises, most likely, 
due to Synchrotron losses at the energy where the acceleration time equals to the energy loss time 
( |de Jager|1996l ): 

7Max. = 4 X 10V5-^/^ (4) 



where / is a numerical constant of order unity which encompasses the details of the acceleration 
process and B is measured in Gauss. In a shock-acceleration scenario this depends on the amount 
of time the particle spends in the downstream and upstream regions ( Piran fc Nakar|20T0 ; Barniol 
Duran &: Kumar [201 Ij ). The total internal energy of the electrons is then (primes denote the co- 
moving frame): 



p-2 



-IrnNerrieC 



e E' 



(5) 



where E'^^^ is the energy of the flow dissipated into internal energy and the ratio of the total energy 
in these relativistic electrons to the total internal energy is denoted by eg. Notice that does not 
necessarily reflect the ratio of instantaneous energy of the relativistic electrons to the total energy, 
which could be much smaller. 
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The magnetic field, B, is assumed to be constant over the entire emitting region. Its energy is: 

E'^ = ^AnR'^ ^ eBE[^, = eE'^, (6) 

where kR/T is the thickness of the shell in the co-moving frame and is the fraction of magnetic 
to internal energy. In this case, for a Baryoninc dominated flow, cb + Ce ^ 1. Whereas, in a 
Poyinting flux dominated jet, one can imagine a situation in which at first the energy is stored 
within the magnetic fields, and only later it is dissipated and transferred to the electrons. Within 
the definitions we use here, this means that eg + < 2. 

Overall, we find that six independent parameters define the synchrotron solution]^ These are: 
the co-moving magnetic field strength, B, the number of relativistic electrons, Ne, the ratio be- 
tween the magnetic energy and the kinetic energy of the electrons, e = e^/ee, the bulk Lorentz 
factor of the jet with respect to the GRB host galaxy, F, the minimum electrons' Lorentz factor in 
the source frame, 7^, and the ratio between the shell crossing time and the angular timescale, k. 
The first three observables: the spectral peak, the peak luminosity and the pulse duration, described 
in ^2.1 reduce the number of free parameters from six to three, which we choose as: (cB/ce, F, A;). 
The other constraints limit the allowed regions within this parameter space. 



3. Detailed method and results 



3.1. The synchrotron equations 



The typical synchrotron frequency and power produced by an electron with a Lorentz factor 
7e are: ( [Rybicki and Lightman|1979||Wijers and Galama|1999[ ) 

2 QeB 



^syn 



27rmf.c' 



Pile 



4 ^ 2^' 



(7) 



(8) 



where F is the Lorentz factor of the bulk and ax is the Thomson cross-section. When the electrons 
have a power-law distribution of energies above some 7^ there are two typical breaks in the spec- 
trum. The first is the synchrotron frequency of a typical electron, z/^ = ^synilm)- The second is 
the cooling frequency, the frequency at which electrons cooling on the pulse time-scale, radiate: 



a^Bmiil + Yo 



\2 ' 



(9) 



Kumar & McMahon 



(2008 1 set the pulse duration to be equal to the angular time scale. We consider here a more 



general approach in which this is only a lower limit. 
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where Yq is the relative energy loss by IC as compared with synchrotron for 7e = 7^ 



The peak of the synchrotron uF,, is at max[iyc, Um] (Sari et al. 1998). We identify this peak 
with the observed sub-MeV peak. For < i^m, the "typical" electrons are fast cooling and rapidly 
dispose of all their kinetic energy, while for > I'm the electrons are slow cooling and the typical 
electrons do not emit all their energy which is eventually lost by adiabatic cooling. We consider 
fast cooling solutions (or marginally fast) for the prompt emission. 



(10) 



The maximal spectral flux, at min[vc, Vm] is given by combining Eqns. [7|[8| 



For z/c < z/m, this is the flux at and the peak of F^. However, in the fast cooling regime the flux 
between Uf. and Urn varies as Fy oc z/^^/^ and the peak of uF^, is at z/^: 

c^arTBxinsNe ,,,, 

^ u,m ^ u,raax\ ) o /( j2 /1 i ^ ' V^^/ 

z/c 3qe4ndl{l + Yo) 

where the factor Xrns/ (1 + ^0) = (t'c/'^m)^^^ accounts for the fact that only Xins of the relativistic 
electrons are instantaneously emitting and only 1/(1 + Yq) of their energy loss is by synchrotron 
(the rest is by the accompanying SSC). In the fast cooling regime, the electrons cool on a shorter 
timescale than the pulse duration and several consecutive cooling processes may be needed to 



account for the comparatively long pulse duration. Eq. 11 can be understood in terms of the 
cooling time of the electrons, tc, the time it takes the electrons to radiate (by synchrotron and IC) 
their kinetic energy: 

We) = . . (12) 

tc(7e) is declining as a function of 7e. Therefore it is maximal for 7^, and after a time tc(7m) all of 
the electrons have cooled significantly. Combining Eqs. 7|[9]and 12 gives: tc(7m)/ip = ('^c/'^m)^^^- 
It follows that: Xj„s/(1 + Yq) = tc{'ym) / tp, which tel 



s us that indeed only the instantaneously 



synchrotron emitting electrons contribute to the peak of uF,^. 

Given the spectrum, the source must be optically thin for scatterings: 

^ Ntotax ^ 

where Ntot = + 2Npairs is the total number of electrons and positrons in the flow. Ntot is 
composed of two components. The first, A'^e, is calculated from the peak flux. Notice that in 
principle there could be additional non-relativistic electrons in the initial flow that although they 
don't contribute to the radiation, they can be important due to their effect on the optical depth. 



We return to this point in {3.2 The second, 2Npairs, is the electron-positron pairs created by 
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annihilation of the high energy photons Nj 
pair creation ( [Lithwick fc Sari| ( |2001[ )): 



pairs 



is derived by considering the optical depth for 



(14) 



where n^u,an is the number of photons in the source with enough energy to annihilate the photon 
with frequency, v. The frequency I'thickO^), is defined as the (source frame) photon frequency at 
which T^^^{h'thick) = 1- Therefore, Npairs equals the number of photons above ythickiX)- 

A more detailed calculation of the pair creation opacity has to take care of the details of the 
radiation field in the source ( Granot et al.||2008[ [Hascoet et al.||2012[ ). There could in principle be 
a numerica l factor in Eq. [T4]of order ~ 0.01 - 0.1 ( [Granot et al.|2008t [Hascoet et al.|2012||Vurm 
et al.||2012| ) This in turn, would lower the minimum requirement on T. However: F oc re, which 
means that the change in the lower limit of V in the exact calculation, will only be of order 0.5-0.7. 



3.2. A synchrotron solution for the MeV peak 



We solve the synchrotron Eqns. described in ^:3.1 and reduce the 6D parameter space de 



scribed in [2.2 to a 3D parameter space defined by (e, F, k). We then consider three constraints 
that limit the allowed regions within this parameter space. In ^|5] we discus further constraints on 
this parameter space that arise from pair creation opacity hmits. 



Solving Eqns. 



and 1 1 we find: 



B = 2x + Y,f/'V-^\ef"k-''\k + lf'^^l%,t-l,F]'l,d2s Gauss, 

X... = 1.1 X 10-^(1 + Y,ry^Ti{er'i^k'/\k + l)-^/v;3^/o^t;2F^,Y.t4^ 

iV. = 1.8 X 10^^(1 + Fo)^/^r-(e)V^A:-V4(^ + l)^/^4oo^;Si^-S.24^ 
R = 3x lO^^Tlik + l)~%,.3 cm, 
7™ = 3600(1 + Yo)-'/'r,{e)-'/'k'/\k + 1)-'''i^I%/J,IfzI^^^^^^^ 

^Ma. = 2.8 X 10^(1 + ro)-^/^rf (e)-vvA:V4(fc + i)-^/V;3l%3F-i/.t4-8 

where F2 = F/lOO, z/p,300 = /ii^peafc/SOOKeV, F_26.2 = i^^,m/(1.5 x lO-^Sgrg sec-^Hz-^cm-^ 
tp^,3 = tp/0.5sec and = di/lO'^^cm. 



(15) 

(16) 
(17) 
(18) 
(19) 
(20) 



^The pairs are created with energies of order half that of the high energy photon. This causes the pairs to have a 
steep power law of random Lorentz factors extending up to very high energies. These pairs radiate by synchrotron and 
contribute mainly at the optical band. Their radiation in the sub-MeV band is sub-dominant. 
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We proceed to consider the limits on the parameter space (e, F, k) that arise from (i) efficiency 
(fast cooling), (ii) the size of the emitting radius and (iii) the opacity. 

Efficiency and fast cooling require Vc < i^m ( §3.1| ), which yields: 

r, < 17(1 + YoY/'%ef''k-'/'%k + l)^/^^4'oW.f ^-26'24'- (21) 



The emitting radius must be smaller than Rdec, the deceleration radius CPanaitescu and Kumar 



2004[ [Kumar fc McMahon|2008HZou et al.||2009i ). The latter depends on the circum-burst density 



profile: a wind with n = {A/iripjr ^ora constant ISM with n Uq. The hmit on the emitting 
region can be written, in these two cases, as: 



R < Rdec 



2.6 X lO^^Egf no ^/^Fg ^/Vax(l, O.Olnl^'^T^^^ E^^^^r2^^^)cm, constant ISM 
1.8 X W^E^^A-^T^'^msixil, 0.36Al^^T^^'^E^l^^Tl)cm, wind environment 

(22) 

where -E54 = Eiso/lO^^erg is the total (isotropic) energy of the burst (and not just the energy 
emitted in one pulse), T2 = T/lOOsec is the total duration of the burst, no is the particle density in 

cm"^ and = A/{5x lO^^g/cm). 

The flow is optically thin for scattering. Typically, we find that Ntot is dominated by the 
pairs created within the flow by the annihilation of the high energy part of the spectrum. The 
main contribution to the number of photons available for pair creation is from IC scatterings of 
synchrotron photons to energies above the pair creation threshold. Solving Eqns. [13 14 and 19 we 
obtain: 

r = 8.7(1 + Yof^T^^er^'k'^^^k + l)^^/«4^ooC/'^-2£4f • (23) 

The criterion r < 1 leads to a lower limit on F. It should be noted that in many models one may 
expect a large amount of passive (i.e. non relativistic) electrons in the flow ( Bosnjak et al.||2009 ). 



In this case the optical depth may become dominated by these electrons. We discuss this option in 

So far we discussed only the synchrotron signal itself, and we haven't yet addressed the SSC 
contribution. We now turn to describe the SSC emission, and explore its effect on the solutions. 



3.3. The SSC component 



Synchrotron is essentially accompanied by SSC. We describe here a single-zone analytic so- 
lution of this component which is expected to reproduce the general characteristics of the true 
solution. However, we note that a full description of the up-scattered flux can only be done with 
more detailed calculations ( Nakar etal.|2009 , Bosnjak et al. 2009 1 and those might differ somewhat 
from the simplified analysis. If a SSC signal is detected, we can measure both the peak frequency 
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and the flux of the SSC component. In this case, we obtain two additional equations that allow 
us to determine the parameters up to a single free parameter which we choose as k. This may be 
the situation in the future, when CTA begins observing in the very high energy range. We return 
to this intriguing possibility in ^ But even if the SSC signal is not detected (as is the current 
observational situation) we can use the non-detection as a limit on the parameter space. Note that 
the signal could be undetected either because it is very weak or because it peaks at sufficiently high 
frequencies and cannot be observed by current detectors. 



The typical SSC photon frequency is given by: 



J^SSC = lm^m- 



(24) 



The IC process is suppressed by the Klein-Nishina effect and the photons are up-scattered only to 
vkn, if the up-scattered photon's energy in the electron's rest frame exceeds the electron rest mass 
energy: 

huKN = T'j^m^c'^ < 'jl.hum. (25) 

This condition is satisfied if: 

e < 1.8 X lO'k'\l + k)'ul,,,tl,F:^,,,d^i. (26) 

Thus, for typical parameters the up-scattered photons are most likely in the KN regime. We intro- 
duce the parameter ( Ando et al.|2008[ ) to take the KN effect into account: 



c 



KN 



' ^SSC ' 



' JmhVTn 



0.22(1 + ro)^(e)iA;4(l + A:)lz/;|ootp;|Fi6.24, if 7^ > 1 



VSSC 

else. 



(27) 



The total SSC flux (integrated over frequencies), Fssc, is related to the total flux ascosiated 
with the synchrotron peak, ~ i'mFu,m, by: 

Fssc 



u F 



= Yo = YC 



KN, 



where: 



Y 



if Ce < Cfi 



(28) 



(29) 



is the Compton parameter ( |Sari et al.|1996| ) . 
This yields: 



0.22e-8/.-8(l + A;)8z/;|oot;|F_WI^^ if 6> l(Fo <>^< 1) 
Yo=i022e-lk~Hl + k)lu;l,t-iFi,Jl,, iil>e >1.7xlO-^k-l{l + k)u;l^^Fi,^,4s(Yo<l<Y) 



_0.17e-U-^(l + A;)7zyJooV3^-Wl8, else (KYoKY) 



(30) 
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As mentioned above, a possible future detection of an SSC component above the LAT range 
can, in principle, measure: ukn and Fssc- This provides two extra equations, allowing us to solve 
for r and e. Using Eqns. 28|29 



we can write: 



1.7 X io-2y^/c4(i+/c)z/;;iot;k' 
1.7 ^ iQ-^Y,Xk'Hi+k)v;l,t;lF 



2 



-26. 2^285 
1 2 
3 ^3 
-26.2"28) 



and 



14 y. 7 p7 



T- 



),300''p,.3 

1 1 



mY-l,k-^{l + k)^v-^^ 

mY,Xk-l{i+k)lv;l,t;%F\^ 



2 

:"28 
1 

:"28 
1 



hVKN ^ 


1 


100Gel/> 






1 


imGeV ' 




hUKN ^ 


1 


IQOGeV ' 





iiYQ<Y<l, 

if Fo (31) 

else (i<ro<>'); 
if ro<>"<i, 

ifY^KKY, (32) 
else (l<yo<>')- 



In addition, detection of a SSC signal above tens of GeV, implies that the spectral breaks 
observed by the LAT are not due to an optical depth. If the particles are shock accelerated then 
these breaks would be associated with the maximal synchrotron frequency: ysyn,max 



4.5 X 



10 r2Hz( |de Jager||1996| ) where / depends on the acceleration mechanism, and is expected to 



be of order unity (Piran & Nakar 2010 [Barniol Duran Kumar 2011 1. If we can also observe 



z/, 



syn,m,ax 



we obtain another equation that would allow us to solve for k and either obtain the full 



solution or rule it out if no solution is found: 

r = ioof-^{hi^, 



syn,max 



/IGeV) 



n f-2\^-l ( hVsyn.max \ ( hUKI\ , 



4.5 X 10-^ V-1 



4 X 



IGeV 



-2 /hvs 



IGeV 



V WQGeV ' 

' hUKN ' 
. lOOGey ' 



p,300' 



'^p,300' 



if ro<>'<i 

HYqKKY 
else {KYqKY), 



(33) 



(34) 



2.2/ '^^o!-l(^^^T^F^)KiMGfy) ''^p,300^Pv3-^-26.2'^28^5 if ^0 < '^^ < 1 

A;=<J 4r%_i(^^^^^)3(^)-fz/,,3o^^^^^ iiY,<l<Y (35) 



IGeV I UOOGeV' 

4 7^-6y^_^(/^^f^)3(_j|^ else {\<Y^<Y). 

As opposed to similar observations in Blazars, current observations in the GeV range do 
not detect a significant high energy component in typical GRBs. Current detections of a GeV 
component are at the level of 0.03 of the MeV emission and upper limits for most other bursts are 
at the level of 0.1. In a few cases (e.g. GRB 090926A ( [Ackermann et al.|20TT] )) an additional high 



energy component has been detected, and may be associated with SSC, however, even in this case 
one does not see a clear GeV peak, like the one observed in Blazars. In this case, the upper limits 
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on the SSC signal further constrain the possible parameter space. First, the photon with largest 
observed frequency, Umax niust satisfy: Umax ^ ^mck, where the latter is the frequency at which 
the flow becomes optically thick to pair creation. Above Umax, the photons could (but do not have 
to) be optically thick for pair creation. The condition Tj^j{i'max) ^ 1 leads to lower limits on T. 
However, this condition is somewhat model dependent ( Granot et al.||2008 , Zou et al.|2009 , 2011 



Hascoet et al. 2012), and the limits on F from this consideration may be alleviated in case of a 
two zone model. We therefore separate the discussion of this limit from the main analysis, and 
return to it later on in ^|5| Second, Beniamini et al. (2011 1 have shown that for a typical GBM 
burst the total flux which is observed in the LAT band (30 MeV-300 GeV) is at most 0.13 of the 
GBM (8KeV-40MeV) flux of the same burst (see also |Ando et al.|2008i[Guetta et al.|2"0TT] ). More 
recent studies, using extra noise cuts applied by the LAT team ( Ackermann et al.|2012 ) give more 
constraining upper limits which are lower by a factor ~ 3. We define tjlat as the upper limit on 
the fractional LAT flux: FLAr/Fy^m ^ 0.1?7lat-i (where tilat-i = //lat/O.I). Below Umax, the 
photons are necessarily optically thin for pair creation and therefore in order for the LAT signal 
to be sufficiently low, either the total up-scattered flux is low (less than r]LAT of the GBM flux), 
or most of the up-scattered flux is at u > Umax and only a small fraction can be observed (for an 
illustration of these possibilities see Fig. [2]). This leads to lower limits on F and on e. 

Let a denote the photon spectral index below Up^ak, i-C- oc u°'. Most GRBs are phe- 
nomenologically well fit by a Band function, and for these a ^ —1. Defining (yy as the fraction 
of up-scattered flux observed in the LAT window, relative to the total up-scattered flux, the ratio 
between the expected flux up to Umax and the sub-MeV flux, is: 



u, 



LAT 



(36) 



IfuxN > ^max and if uF„ is rising below ukn (as happens for a ^ — 1) then the total flux observed 




Fig. 2. — Two schematic possibilities allowing the SSC component to be compatible with LAT observa- 
tions. 
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in the LAT window is dominated by Umax- Cw can be written as ( Ando et al.|2008 1 



Cvy ~ < • r 1 (37) 

^max ' 

Plugging in ukn (i-e. assuming ukn < ^ssc) yields: 



(38) 



where hVrnaxA 



IGeV 



We require that the SSC contribution will be sufficiently low to agree with the LAT observa- 
tions. We separate the possible solutions to two cases. 

If e > 2.4A;-V7(i + A;)3/7^^-5A^-f = Yo < O.It^lat^^ and the 

total up-scattered flux (and not just the fraction within the considered band) is below the observed 
limit. Otherwise, the total ratio of up-scattered to synchrotron flux is larger than 0.1r]LAT-i, but it 
peaks at u^n ^ t^max and thus may be partially absorbed. As T increases, the up-scattered flux 
peaks at higher frequencies. This means, that for large enough values of F, the up-scattered tail 
below Umax is small enough that it is compatible with the observations. Therefore, in this regime 
of e, one can obtain a lower limit on T in terms of r]LAT arising from the requirement that the total 
up-scattered flux up to umax be less than the flux observed by LAT: 

2.25a-4.27 1/9 2a + l 2a + 5 . . 6a + 15 - JS'^'^A , - m?^I^a\ Jn^^A^ 

> io-^^r^r/T-i^ilie«fc"^(i+^)^^p;or^^ (39) 

Notice that these are conservative limits as we only use the spectral range below Vmax where we 
know the conditions are optically thin. One must bear in mind, that efficiency considerations limit 



the amount of energy that may be carried by the up-scattered flux by virtue of Eq. 41 This limits 
e but not F. 



3.4. Results 

Figs. [3| |4] depict the different GRB parameters superimposed on the allowed space in the 
(e, F) plane arising from the above conditions. We plot here results for A; = 1 and k = 10. The 
condition that the SSC flux resides below the observational limits within the frequency range up to 
the maximal observed photon energy is plotted for Umax = 71GeV (source frame) corresponding 
to the highest (prompt) energy photon to date, from 0809 16C. They are drawn for two cases of the 
lower spectral slope: a = —1 and a = —1.5. We choose to take the most extreme values for this 
limit, in order to show that even in this case, there is still reasonable parmeter space in the (e, F) 
plane. We return to this extreme case in greater detail in ^ 
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Several characteristic features can be seen in these figures. First, R is directly related to the 
pulse duration and bulk Lorentz factor via Eq. |2[ and it is independent of e. For the canonical 
observed parameters we use here, R spans two orders of magnitude, 10^^A;^^/^(A; + .3cm < 

R < 10^"^ k^^^^{k + l)^/'^tp 3cm, and it is relatively large. B is the most fluctuating parame- 
ter spanning almost 4 orders of magnitude: 3k^^/^^{k + l)^^^^t'p,3ooVr3'^-26.2'^28 Gauss < B < 
10'^k^/^^{k + 300V 3-^-26 2*^28 Gauss. The lines of constant B are almost parallel to the 

= line, which means that the value of B is almost a direct representative of x^s, the 
fraction of instantaneously emitting electrons relative to the total number of electrons. The al- 
lowed range for 7^ is: 3 x 10^k^/^%k + l)"^^/^®z^p%tJ5Fr2/.t^^28^^^ < Im < lO^k^^^^k + 
]^-)-3/i6^^/4^^i/2^-V4^-^i/2_ jj^ggg ]^jgj^ values of 7m (compared with nip/me) are extremely sig- 
nificant for GRB models in which the electrons are heated by shocks and the initial energy resides 
in protons. In these cases, it is necessary that only a small fraction of electrons will be heated to 
relativistic velocities in order to allow them to reach such high energies (Daigne & Mochkovitch 



1998[ Bosnjak et al.|2009 ). We show this explicitly for the internal shocks model in ^ In addition, 
we observe that i^kn is expected to lie in the range: 250k^^^{k + l)'^^^'^'^p,3oo^pfl-^r26Trf28^^^GeV < 
ukn < l-7k~^^^{k + l)'^^^'^p 3oo^p'^3-^r2Y2'^28^'^^TeV. Therefore, even with no knowledge of the 
LAT observations, we could have expected that the SSC does not peak in the GeV but typically at 
least two orders of magnitude above. This is a direct consequence of the large values of 7^ required 
for this solution. A^^, is expected to lie between: 105°/c-i/S(fc + l)-^/24^3/4^^y2^5/4 ^^5/2 < < 

10^^fcV8(fc + l)-3/8,3/^4^,V2^ 

Even for the highest observed photon energy of 71 GeV, we see that the lower limits on T that 
arise from the SSC flux in the GeV are less constraining than the lower limits that arise from the 
optical depth. However, the SSC limits become very strong for small values of the lower spectral 
slope, as expected in the fast cooling synchrotron regime. For the expected a = —1.5 (which is 
found for 10% of the GRBs in the GEM and BATSE samples) e > 3 x 10"^ for = 1 (e > 2 x 10"^ 
for A; = 10). In addition, very negative values of a, push solutions with relatively low e towards 



the marginally fast solution {Uc ~ i^m) as discussed in §3.6[ By Eq. |23| the lower limit on F varies 
as k~^^^^{l + i^) 25/48 therefore increases for A; ^ 1. The upper limit that arises from < Um 
scales as ^-^/^^(l + K)^/^^, and also increases for A; ^ 1 . similarly the allowed parameter space 
mcreases withQA; § 1. 

The main effect of increasing k is to increase the allowed values for F. The allowed range 
for F is roughly 150 < F < 3000 for A; = 1 and 300 < F < 6000 for k = 10. This range 
is narrowed down by considering the pair creation opacity limits and it will be addressed again 
in ^|5| Interestingly, increasing k does not change significantly other parameters of the model. 



"*The less rigid upper limit on the parameter space arising from R < Rdec also increases with k, either as (k + 1)'^/^ 
for an ISM or as (fc + 1)^/* for a wind environment. 
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R, B, andNe- R oc T'^/k implies that the lower limit on the radius scales as k^^/^{k + 1)^/^'* 
while the upper limit scales as k~^/^{k + 1)^/^. Thus, even for relatively large k one expects a 
comparable range of emission radii to the one obtained for A; = 1. For the magnetic field, the lower 
limit scales as k^^^^^(k + 1)^/^^ while the upper limit depends on k^/^^{k + 1)"^/^^. Again, the 
allowed range of B varies very weakly with k, and 3 Gauss and 10"^ Gauss are strict lower and 
upper limits. The range of electrons Lorentz factors varies with a lower limit oc k^^^^{k + i)-ii/48^ 
and an upper limit oc k^^^^{k + 1)^^/^^ yielding a strict upper limit of ~ 10^ (obtained for k = 1) 
and a weakly varying lower limit of order 3000. Finally, the lower limit on the number of electrons 
is oc k^^/^ik + 1)^^/^^, while the upper limit is oc k^^^{k + implying a strong upper limit 

on the number of relativistic electrons ~ 10^^. 



3.5. Energetics and efficiency 

We denote by Eint the fraction of the overall energy in the flow, Etou that is dissipated into 
internal energy. In turn, part of this energy, Egyn, is radiated by synchrotron to produce the sub- 
MeV peak. Thus, we define: 

Egyn = VsynEint = VsynVintEtot = VtotEtot- (40) 

rjint is the conversion efficiency of bulk to internal energy in the prompt phase, i.e. the efficiency 
of the energy dissipation process, and rjsyn < 1 is the conversion efficiency from internal energy 
to radiation. GRBs must be highly efficient and must emit a significant fraction of the total ki- 
netic energy in the sub-MeV to avoid an "energy crisis", i.e. rjtot cannot be much less than unity 
dPanaitescu and Kumar|2002HGranot et al.|2006HFan et al.|2006D . 



The value of rjint is determined by the energy dissipation process and is therefore not the 
subject of this paper. We remark, however, that in many models 'qmt is expected to be relatively 



low. For example, in the internal shocks model it is expected that r]int < 0.2 ( jKobayashi et al. 



l997||Daigne fc Mochkov"itchl[T998| pelborodov|[2000t [Guetta efaLjllOOT ]). Notice, however, that 



an exception is possible,in case the same energy can be reprocessed to internal energy more than 
once leading to r]int > 1 ( Kobayashi et al.|[2001 ). It is unknown how natural this scenario is and 



we leave the discussion on this possibility for future studies. We conclude that quite generally rjsyn 
must be smaller than unity. 

Two efficiency criteria determine the efficiency of the synchrotron signal. First, the energy 
radiated by synchrotron is limited by the energy of the relativistic electrons. Therefore, we focus 
here on fast cooling to avoid an extra loss in the efficiency. This means that r]syn ^ and therefore 
eg cannot be too small. Accounting for SSC losses in addition to synchrotron losses, this translates 
to an upper limit on e. On the other hand, e is also bounded from below by the requirement that 
SSC loses are not too large. = Yq is a function of e, independent of the bulk Lorentz factor. 
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Fig. 3. — The range of GRB parameters for typical observables (tp.3,F_26.2>c^28>2^p,30o) in the (e, T) plane 
for k = 1. Plotted here (from top left in clockwise order) are: i?, B, jm, Ne- Colours depict the values of 
these parameters. The conditions r < 1 (bottom horizontal line) and Uc < Vm (top curved line) impose 
strict limits on the parameter space. The area between them is available for Synchrotron solutions. The 
requirement that the SSC signal is below the observational limits in the GeV range, leads to lower limits on 
r (e) which are depicted by the blue circles for the typical lower spectral index (a = — 1) and by green X's for 
the expected slope in the fast cooling regime (a = —1.5). The sharp cut-off in these limits corresponds to the 
transition between ^^^^ =Yn< 0.1?7LyiT-i where the total up-scattered flux is below the observational 
limit, to ^-^^^ = Yn > O.lriLAT-i where the up-scattered signal should peak at high enough energies so 
that the flux residing in the observed band is sufficiently low to account for the observations. These are cases 
1 and 2 in Fig. [2] correspondingly. Two other conditions R < i?rfec(Wind) and R < Rdeci^SM) impose 
softer limits (the areas below these lines are allowed for each case) on the possible solutions as they both 
depend on additional unknown parameters (A^, and no). 
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10^ ^GeVlimit(a=-1) 
■ GeVlimit(a=-1.5) 
-H.n„(ISM) 



Fig. 4. — Same as Fig. |3]but for k = 10. 



10^ —GeV limit (ii=-l) 
■ GeVlimit(a=-1.5) 



(o)l-og(vJ 



-- T = 1 

lO^ -^GeV limit (u=-1 
■ GeV limit (a=-1 
-R=R^^^(ISM) 
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but in the KN regime it may still depend on k and on the various observed parameters. As seen 



from Eq. 30, apart for the peak frequency and pulse duration (which can span over a large range 
of values), all of these dependencies are very weak and have no significant effect on the overall 
efficiency. We examine the dependence of r]syn on Up and tp in Fig. [5} 

Combining all these efficiency factors, we get the overall efficiency: 



(41) 



As rjint alone is already expected to be relatively low, this limits rfsyn ^ 0.1 in order to assure 
that the overall efficiency satisfies: rjtot > 10^^. Thus, we must have: 10^"^ < e < 10 (see Fig. 
[5]). Notice that this result is independent of the bulk Lorentz factor, and it is weakly dependent on 
k and the observed parameters. Out of these, the strongest dependence is on Up. Vp ^ lOOKeV 
(z/p ^ lOOOKeV) increases (decreases) the lower limit on e to about 10"^ (10^^). Although low 
t'peafe results in low efficiencies, such bursts are usually observed with lower luminosities, and the 
energetic requirements for such bursts are not unreasonable. We return to this point in more detail 



in ^: 3.7 The efficiency depends more weakly on the pulse duration, but the latter spans at least three 
orders of magnitude and thus can cause noticeable change. Taking tp ^ 5 x 10~^sec {tp ^ 5sec) 
yields e > 3 x 10~^ (e > 10"^). Later on, we find that high magnetization is favored by other 
reasons as well. 
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Fig. 5. — An upper limit on the efficiency when the sub-MeV peak is ascosiated with fast cooling syn- 
chrotron emission, as a function of e. The electrons must carry a significant amount of the energy, but too 
strong an amount will lead to a large amount of energy released at frequencies far above the peak by the 
KN process. The left panel shows the dependence on Up while the right panel shows the dependence on the 
pulse duration. These plots depict upper limits on the efficiency, as they do not include the efficiency of the 
energy dissipation process or energy losses in components other than magnetic fields or electrons. 



3.6. The spectral shape 



The solutions given in 3.2 consider only the peak frequency and the peak flux, disregarding 
the spectral shape. In the fast cooling scenario (considered here) the slope of the photons' number 



spectrum above the peak of z/Fj, is (z/'^ oc i^^^P+'^y^) ( |Sari and Piran 1997a). This is consistent 
with the average observed slope (/3 = —2.3) for p ^ 2.6 which is approximately the expected 
value for p in case the electrons are accelerated by the Fermi mechanism ( Achterberg et al.||2001| ) 
(Although it is not clear that the spread of (3 is consistent with the expected spread in p (Shenj 
2006 [ |Starling||2008[ |Curran|[2010[ )). However, the synchrotron fast cooling slope below the peak 
of lyF^ iy" oc i/"^-^) ( [Cohen et al.1|l997i |Sari et al.|[T998l [Ghisellini et"aLl|2000K is inconsistent 
with the average observed slope a = —1 and 90% of bursts have steeper slopes than this ([Preece 



et al.|2000| [Ghirlanda et al.|2002| |Kaneko|2006| |Nava etal.|2011D . Furthermore, the lower energy 
slope in about 40% of GRBs is a > —2/3 ( Nava et al.|2011 1 which is impossible in a synchrotron 



scenario, even in the case of slow cooling. This is known as the synchrotron "line of death" (Preece 



et al. 1998 1. An extra concern, is that the break observed at fpeak is often too sharp to be compatible 
with the smooth transition expected from the synchrotron model ( Pelaez et al.|1994| ). 



The synchrotron "fast cooling - line of death" problem, the more serious of the two, can be 
partially resolved if the electrons are marginally fast cooling ( Daigne et aLpOlT ) with i/^ ^ u^- 
This requires, of course, an additional fine tuning mechanism that keeps this condition satisfied. To 
examine the possible parameter phase space with a "marginal fast cooling" we define a marginally 
fast solution as one that obeys one of two criteria. Either it is a fast cooling solution (z/^ > fc) 
that obeys: < St'c (If were much larger, there would be a large range of frequencies where 
a = —1.5 and the "fast cooling - line of death" problem would not be solved); or it is a a slow 
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cooling solution (Uc > i^m) (where a = —2/3 is always valid so long as we associate u^, and 
not Uc, with the observed sub-MeV peak) which is sufficiently efficient, with: r]syn > 0.1. This 
family of solutions, is seen in Fig. |6]and it occupies a significant area in the parameter space. The 
break in the slopes of the parameters in the (F, e) plane is due to the switch from the fast to slow 
cooling regime. The upper curve is the limit on the marginally fast solution due to efficiency. These 
lines cross the = J^m line at e ^ 10^"^ and e ~ 10 which is where the fast cooling efficiency 
falls below 0.1 (see Fig. [5]). The marginally-fast solutions are characterized by tc(7m) ~ tp, and 
therefore Xmt ~ 1 and the instantaneous number of emitting electrons is about the same as the 
total number of relativistic electrons. The allowed region has high values of F and 7^: F > 700, 
7m ~ 10^, relatively small values of N^,: ^ 10^°, a weak (and almost constant) magnetic field: 
5^10 Gauss and a large emission radius i? > 3 x 10 
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cm. 



A canonical wind solution is ruled 
out for these solutions. Another partial solution to the "line of death" problem, can be achieved 



in case SSC occurs in the KN regime, as was shown in ^3.3 to be the generic case. In such cases. 



detailed modeling, shows that a can be increased from -1.5 to -1 ( |Derishev]|2003[ |Nakar et aL 
2009|[mng|2009|[Daigne et al.|20lT||Bamiol Duran et al.|2"0T2l ). However, this solution requires 



very low values of e and with e > 10 as required by efficiency considerations (see ^ 3.5 1, one is 
still limited to a < -1.2. 




Fig. 6. — The range of GRB parameters for a marginally fast solution in the (e, F) plane for = 1, ^ = 1. 
Plotted here (from top left in clockwise order) are: R, B, jm, N^. Below the dashed line {vm = 3fc) there is 
a sufficiently large frequency range between and leading to a = —1.5 below the sub-MeV peak which 
is incompatible with the marginally fast solution. Above the black dot-dashed line, the radiative efficiency 
falls below 0.1. In addition are plotted the lines of r]syn = 0.5, 0.3 in grey dot-dashed lines. 
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A complete solution to the "line of death" problem within synchrotron emission can be 
achieved if the synchrotron self absorption frequency, vssa, is sufficiently close to Vm- Extrap- 
olating the fast cooling synchrotron flux below and equating it to the black body limit ( |Granot 
et al. 12000] ) we can estimate ussa'- 



2u 



SSA 



'l{l'SSA)^me& 



(42) 



where '^{yssA) is the typical thermal Lorentz factor of electrons radiating at vssa- Eq. 42 yields a 
(source frame) self absorption frequency of: 



Yssa_ 



9.2 X 10-5(1 + ro)'/''r-'/=^(e)Vi2A;-i/i2(A; + l)ii/i2^J/4^-f 



Typically, pssa is about 4 orders of magnitude bellow Up and it depends weakly on the various 
parameters. It is therefore, unlikely that self absorption has any significant effect on a. 

A final possibility that should be mentioned in the "line of death" context is that a may be 
affected by the presence of a sub-dominant thermal component which is expected to occur in the 
fireball model at a level of a few percent of the total prompt energy ( [Meszaros fc Rees|2000[|Daigne 
fc Mochkovitch|[2002t |Nakar et al.||2005[ ). Detailed spectral modeling, shows that adding a weak 
sub-dominant thermal component usually causes a to decrease ([Guiriec et al.|2011j|2012[). 



3.7. The Distribution of Epeak\Lpeak 



There have been many claims that the distribution of Epeak is rather narrow (Band 1993 , Mal- 
lozzi et aL||1995[ [Brainerd et al.||1999| |Schaefer||2003[) . Typically Ep^ak values are in the range 



150 KeV < Epeak < 700 KeV (Kaneko 2006). It is not obvious that this observation reflects 



the intrinsic properties of GRBs and is not caused by some selection effect (Shahmoradi 2010[ ). 
Softer bursts have been detected by Beppo-SAX and HETE2, and Fermi is finding much higher 
Epeak bursts. Howcvcr, if true, this result is surprising in the general framework of the Synchrotron 
model, where one expects Ep^ak = oc T'j'^B and it is not obvious why this specific com- 
bination of parameters will remain fairly constant between different pulses and different bursts. 
Furthermore, there have been claims that Epeak may be correlated with the peak luminosity, Lpeak, 
( [Yonetoku et al.|2010] ) and that therefore the distributions of both are not statistically independent. 

To examine this issue we examine four characteristic bursts in the {Epeak, Lpeak) plane: 

1. Low Epeak (100i^'e\^) and low Lpeak (3 x 10^^ erg /sec) (low end of Yonetou relation). 

2. High Epeak (lOOOKeV) and high Lpeak (1-5 x lO^^erg/sec) (high end of Yonetoku relation). 
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3. High Epeak (lOOOKeV) and low Lp^ak (3 x 10^^ erg /sec) ("above" the Yonetoku relation). 
None of the bursts with known redshift belong to this category, but this may be an observa- 
tional effect dNakar and Piran|2005i |Band|2005D . 



4. Low Epeak (lOOKeV) and high Lp^ak (1-5 x lO^^erg/sec) ("below" the Yonetoku relation). 
Such bursts would be relatively easier to detect (many photons) but are not seen in observa- 
tions. It is therefore believed that they do not exist. 

We are interested in the range of parameters ascosiated with synchrotron solutions for these bursts. 
We note that Lpeak oc F^^^^^UpeakCp , so that there is generally more than one way to realize a 
specific value of Lp^ak with the different observed parameters. In the following discussion, we 
choose to take d to be constant and vary F^^^^^ . In terms of the allowed ranges for the synchrotron 
parameters we do not find any significant difference between the four typical bursts (see Fig. |7]for 
an example of 7m in each of these bursts). The main difference is a slight change in the limits of 
the allowed region corresponding to Uc = and r = 1. This can also be seen in Fig. [8] where we 
plot those limits in the (e, fpeak) plane assuming that Upeak and Lp^ak follow the Yonetoku relation. 

In Fig. |9]we plot the radiative efficiency and minimal luminosity requirements of these bursts 
(as described in §3.5| ). A type-4 burst is the most inefficient and has the largest intrinsic power. 



Its synchrotron efficiency peaks at less than 0.4, and it could very likely be less, depending on the 
value of e, while its intrinsic power is at least 4 x IQ^^erg / sec. This may pose an explanation as to 
why practically no soft - high luminosity bursts are seen (i.e. bursts below the Yonetoku relation). 
The energetic demands for such bursts are simply too high. The other class of Yonetoku relation 
violators, are type 3 bursts. For these, efficiency is high and they cannot be ruled out by the same 



argument. However such bursts (hard and weak) have only a small number of photons (Nakar and 



Piran||2005 , Band||20()5 ) and are observationally hard to detect. This might cause fewer bursts of 



this kind to have a detectable redshift, leading to a natural selection effect. 

Unless, typical bursts have very low e, the observed luminosity is much more important than 
^peak in order to determine the intrinsic power. Therefore, it is possible to produce bursts with 
low Lpeafc and Vpeak and we conclude that this model does not naturally predict the claimed narrow 
distribution of Upeak by efficiency arguments alone. As mentioned above, this may actually be con- 
sidered as an advantage of the synchrotron model, as it is quite likely that the true peak distribution 
is not as narrow as previously thought. In this case, the argument would be reversed. Synchrotron 
would be a very natural way to produce such large peak dispersion. 
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Type 1 burst: E^IOOKeV. L-3 lo^^erg/sec 




Type 2 burst: E^=1000Key, L=1.5 lo^^erg/sec 



(o)Log(7j 



-GeV limit (a=-1) 
■ GeVlimit(ra=-1.5) 
..R=R,(1SM) 



"—GeV limit (a=-1) 
■ GeVlimit(a=-1.5) 
-R=Raec 



Typ6 4but^: E =100KeV, L=1.5 lo^^erg/sec 



Fig. 7. — The range of 7m (shown by different colours) for four types of bursts in the (e, F) plane for k = 1, 
^ = 1. The conditions r < 1 (bottom horizontal line) and Uc < t'm (top curved line) impose strict limits on 
the parameter space. The area between them is available for Synchrotron solutions. The requirement that 
the SSC signal is below the observational limits in the GeV range, leads to lower limits on r(e) which are 
depicted by the blue circles for the typical lower spectral index (a = — 1) and by green X's for the expected 
slope in the fast cooling regime (a = —1.5). Two other conditions R < i?rfec(Wind) and R < iJ^ec (ISM) 
impose softer limits (the areas below these lines are allowed for each case) on the possible solutions. 



4. A synchrotron solution within the internal shocks model 

We divert here from the main theme of the paper, in which we consider a "model independent 
synchrotron emission" and examine specific additional constraints that arise when considering a 
situation when the energy is stored in the bulk kinetic energy of the flow. This is the case in the 
common internal shocks model, but is also the case in a Poynting flux dominated model if the 
magnetic energy has been converted to kinetic energy before the radiation process begun. 

The kinetic energy of the initial flow is dominated by protons Denoting the number of 
protons by Np, we have: Np = Ntot {'^Ne+ = Ntot for a pair dominated plasma). The total 
(bulk+internal) energy needed to produce a combined shell traveling with a bulk Lorentz factor F, 



^In principle, many pairs could be created after the shocks collide and the internal energy is released. This could 
cause the combined shell to be dominated by pairs, but has no effect on the following consideration which is related 
to the initial composition. 
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10 



^peak 



Fig. 8. — Limits on the allowed parameter space due to Um ^ i^c and r < 1 in the (e, I'peakX) space 
assuming that Vpeak and Lp^ak follow the Yonetoku relation ( Yonetoku et al.|20T0 l. Fast cooling optically 



thin synchrotron solution an arise in the area between the two surfaces, 
in this scenario, is: 



Hot 



gTNtotmpC^{x —for pairs) (44) 



where the factor of g ^ few arises from the fact that the colliding shells must initially have 
different bulk Lorentz factors in order for internal shocks to occur (the actual number depends on 
the amount of variations in the bulk Lorentz factor along the flow). The minimal value of g (the 
most efficient) arises when all the initial energy is transformed into two channels: internal energy 
of the electrons and energy in magnetic fields. In this case g = 1/(1 — r]int)- However, if some of 
the energy remains in thermal energy of the the protons, g would be somewhat larger. 



Eq. 44 yields the energy of the emitting electrons: Ee = egEint = e(,r]intEtot in terms of 
the initial energy stored in the protons (or pairs). At the same time the synchrotron conditions 
determine the number of emitting electrons, N^, and their Lorentz factor, 7^ and this determines, 
in turn, the energy of the emitting electrons via Eq. [5j Comparing the two expressions we find 
that the number of emitting electrons is significantly lower than the total number of electrons. We 
define ^ as the fraction of the synchrotron emitting electrons in the flow: ^ = N^/Ntot- A schematic 



picture of the electron number distribution is given in Fig. 10 
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Fig. 9. — The maximal radiative efficiency of the four different types of bursts. Soft, high luminosity bursts 
are the least efficient and require huge energies at the source. This may explain the lack of such bursts. 



p — 2 

Using Eqns. [T7l[T9l|44l|45l = teEint = teVintEtot becomes: 



(45) 



tot 



p- 12 0.2^ " 



X — for pairs). 

nip 



(46) 

The reason for these low fractions of emitting electrons is clear. The synchrotron solution requires 
relatively large electron Lorentz factors, of order 3 x 10^ < 7^ < 10^. On the other hand if the 
available energy eerjmtEtot is distributed to all the electrons (whose number is dictated by the initial 
number of protons) than the typical Lorentz factor would be of order gee'^mt (p— 2) /(p— 1 ) {nip/ m^) 
(or gteriint{p ~ 2)/(]9 — 1) for a pair-dominated flow). Unless g is unreasonably large, this Lorentz 
factors is too small by several orders of magnitude. To overcome this problem we must conclude, 
following Daigne & Mochkovitch (1998]); Bosnjak et al. (2009) that within the internal shocks 
model only a small fraction of the electrons are accelerated in the shocks. It should be noted, 
that the conditions needed for this solution, are very different than the situation achieved in recent 
PIC simulations. This is because, in order for the solution presented here to work there must be a 
gap between the Lorentz factor of the thermal (i.e. non synchrotron emitting) electrons and 7^ as 



depicted in Fig. 10 The ratio of the number of emitting electrons to the total number is plotted in 
Fig. |4](Fig. [TT]for pairs), for the limiting case of + eg = 1 which provides an upper limit on the 
actual value of ^(F, e^, Ce)- Large values of F lead to ^(F, e^, ee) ^ 1. For F > 100 (as found in 
12|), ^ < 4 X 10-2 (^ 2 X 10-5 for pairs). 



In case ^ is very small, there are many "passive" electrons in the flow which may dominate 
the optical depth. is defined such that for ^ < the passive electrons dominate the total number 
of electrons: 



1.4 X 10-^fl 



ro)^/^(e)^/^fc^/^(l + A:)-^/«4«ooC3^-2Y.2rf28^/^ 



(47) 
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In this regime: 

r = 1.1 X 10-^(1 + Yof'%%ey/'k''^\l + ky'/'r'^^^^^^^^^ (48) 




Fig. 10. — The electron distribution is composed of three ingredients: the power-law distributed electrons, 
emitting synchrotron at the sub-MeV peak, the thermal "passive" electrons which are not accelerated to high 
Lorentz factors and the pairs created by annihilation of high energy photons (these too, are not energetic 
enough to contribute to the sub-MeV spectrum). Depending on ^, either the second or the third components 
dominate the total number. 




Log[^{r, es, Ce)] of a proton dominated flow within the internal shocks scenario The results are plotted for 
g = 2, rjint = 0.2, p = 2.5 and + eg = 1. This is an upper limit on the actual value of Log{£) for real 
values of eB, eg. Even for modest values of F, ^ ^ 1 which corresponds to a small fraction of relativistic 
electrons in the emitting region. 



5. pair creation opacity limit (photons above v^ax) 



The limit r^,^(z/„ 



^ 1 leads to a lower limit on F (Woods and Loeb 1995, Piran 1999 



Lithwick k. Sari|2001 1. Eq. |23] shows that in the typical case, most of the electrons in the flow are 
not emitting the synchrotron signal, but are produced by pair creation. Therefore, the limit on F 
from the requirement T^^'^{h'max) ^ 1 is always more constraining than that given by Eq. 23 In 



order to find r^,^(z/„ 



we calculate n 



>max,an 



, the number of photons in the source with enough 
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10 



10 



10 



Fig. 11. — Same as Fig|4]but for a pair dominated flow, i.e. Log[£,[T, cb, ^e)] for g = 2, rjint = 0.2, p = 2.5 
and es + ee = 1. For reasonable values of F, ^ < 2 x 10~^. 

energy to annihilate the highest frequency photon, Urnax by integrating over the spectrum. This 
yields a lower limit on F: 

1 1 2/3-2 

T, > (10-(^-6+^-«/')|l + /3|z/r.,i(l + fc)-V;,^oVp,.3i^--2V2rf28')^- (49) 

This limit is independent of e. 

These lower limits on F replace those used in Figs. |3| |4] in case the high energy photons 
originate from the same site as the sub-MeV emission. F > 300 is required for typical parameters, 
a factor of two larger than the lower limit imposed by considering the optical depth for scattering 



of synchrotron photons. Fig. 12 depicts the allowed region in the (e, F) plane, arising from the 



lower limit on F discussed here and the upper limits given by Eqns. 21 22 The results are plotted 
for the most constraining GRB to date, GRB 0809 16C, for which humax = 71GeV (source frame) 
was observed. For this burst, we find very high values of the Lorentz factors: 10^ < F < lO'', 
10^ < 7m < 3 X 10^, and a large emitting radius 4 x lO^'^cm < i? < 3 x 10^'^cm that is 
only marginally compatible with the canonical wind parameter. Even for this extreme case, the 
synchrotron solution is possible. In fact, the allowed parameter space may be even larger in case 



the high energy emission originates from a different zone than the sub-MeV radiation (Zou et al. 
20091 [20TT|rHascoet et al.|20T2l ). 



As mentioned earlier, given that no opacity break is seen in the spectrum up to Umax, the 
frequency in which the system becomes optically thick to pair creation, Uthick obeys: I'thick > ^max- 
So far, we still have a large uncertainty in the actual value of ytuck- In what follows, we argue that 
this uncertainty can be reduced. The higher uthick lies within the range [i/max300GeV] (300GeV 
is the upper limit of the LAT detector), the implications of the low LAT detection rate becomes 
more constraining on the SSC peak (it should either peak at higher frequencies or be weak enough 
to account for observations). We therefore plot the combination of this pair creation opacity limit 



and the limit from the SSC flux in the GeV (discussed in ^3.3 1, in Fig. 13 for k = 1 and k = 10. 



These lower limits on F are more constraining for lower values of a (causing a smaller loss by 
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extrapolation of the SSC from its peak to the GeV) or for low values of e for which the SSC flux 
is larger. Quantitatively, this is described by Eq. 39 We see that in case umck lies at the highest 
frequency observed by LAT (or above) and a = —1.5, the SSC limits rule out e < 5 x 10^^ for 
A; = 1 (e < 5 X 10-2 for k = 10). 



6. Implications of CTA observations on the Synchrotron model 

It is interesting to explore the implications of future GRB observations by the CTA, in the tens 
of GeV to tens of TeV energy range, on the synchrotron model. The universe is optically thick to 



radiation above a few hundred GeV (Nikishov 1962 Gould et al. 1966) and it will therefore be 



difficult to observe GRBs from typical redshifts z ^ 1 in this energy range. However, the CTA is 
expected to have better detectability in the tens of GeV range than LAT by at least four orders of 
magnitude C CTA consortium ||2010[ 2011). There are three possible scenarios regarding observa- 



tions of the prompt emission by CTA, schematically drawn in Fig. 14 Least intriguing is a null 
detection. This could happen either because the system is optically thick to intrinsic pair creation at 
lower frequencies than the CTA band (i.e. the breaks in the LAT spectrum are due to optical depth 
for pair creation, and the fact that we do not see anything with CTA is trivially expected) or be- 
cause the breaks in the LAT spectrum are intrinsic spectral breaks due to the maximal synchrotron 
frequency and there is simply a very low SSC signal at the CTA band. The latter case could provide 
significant constraints on the parameters of the model. However, since we cannot rule out a large 
optical depth, we won't obtain new constraints on the model just from a null detection. Second, 
there might be a weak signal in the tens of GeV range which is, however, significant enough to 
show that the system is optically thin at least up to the CTA band. In this scenario we may asso- 
ciate the LAT breaks with the maximal synchrotron frequency and may obtain T according to Eq. 



33 up to the numerical factor / which is expected to be of order unity. In addition, the improved 
sensitivity of CTA in the tens of GeV range (relative to that of LAT), could possibly allow us to rule 
out a wide range of SSC solutions (the SSC component should be both weak and peak at higher 
frequencies than u^ax for the SSC solution to be compatible with such observations). Effectively, 



the improved sensitivity would increase the lower limit on F given by Eq. 39 , by replacing Umax 
by the maximal photon energy observed with the CTA and replacing the fractional flux detected by 
LAT, r]LAT with the corresponding fraction r]cTA that could be lower by up to two or three orders 
of magnitude (depending on the actual flux that would be detected by CTA and on the new Umax)- 
This would narrow down the allowed parameter space and would produce lower limits on F and 
e. For instance, a detection of even a few photons in the tens of GeV range, could lead to e > 0.3 
and F > 500. Finally, a strong signal, with a fluence larger than the extrapolated LAT emission 
(i.e. with approximately l^)'"^ ergs /cm?), detected by the CTA would imply a detection of a new 
component in this energy band. Such a component, if it ascosiated with the same emitting region, 
can naturally be explained as the SSC peak If the detection is strong enough, it would allow us to 
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Fig. 12. — The allowed region in the (e, T) plane for GRB 080916C, with k = 1, ^ = I, hvmax = 71GeV 
(source frame). Plotted are: 7^ (left panel) and R (right panel). Colours depict the strength of these 
parameters. The conditions t^^j < 1 (bottom horizontal line) and Vc < Vm (top curved line) impose strict 
limits on the parameter space. The area between them is available for Synchrotron solutions. Two other 
conditions: R < i?dec(Wind) and R < Rdec(iSM) impose softer limits (the areas below these hnes are 
allowed for each case) on the possible solutions. 




Fig. 13. — Lower limits on T from the combined constraints on up-scattered flux and pair creation opacity 
shwon along with the upper limit on F originating from Vc = Vm. for various values of the lower spectral 
slope {a) and the photon energy at which the system becomes optically thick to pair creation {Etuck)- The 
area above between the dashed or dotted lines and the solid line is allowed by observations. Left panel: 
k = 1, right panel: A; = 10. 
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measure both the peak flux and the peak frequency (see ^ 3.3 1. This will allow, in turn, to determine 
ukn and Fssc and eventually to estimate also T and e, and hence all other parameters of the sys- 
tem up to k. In case there values turn out to be outside the limits on the parameter space (z/^ > t^c. 



T > 1 etc.) described in ^3.2 these measurements could possibly rule out the synchrotron solution. 
In addition if we can measure ysyn,max we may obtain another equation and obtain the full solution 
of the model. 





LAT limit ~0.1-1TeV 




LAT limit -0.1-1GeV 


10^ 

CTAJimit ~0>'lTf!V 


-10* 









LAT band CTA band" V 



Fig. 14. — The three possibilities for CTA observations, described in the text. Case A: No signal in the 
CTA - no significant costraints on the parameter space (as this could be due to optical depth). Case B: a 
weak signal at tens of GeV, where the CTA is much more sensitive than LAT, leading to strong constraints 
on SSC. Case C: a detection of a peak (in vFy) in the CTA band, leading to a full solution of the SSC model 
or ruling it out. 



7. Conclusions 

We examined the allowed parameter space for the synchrotron emission for the prompt phase 
of GRBs. We considered individual pulses as the building blocks and demanded that the system is 
fast cooling with the peak {vFy) frequency comparable to z/^ so that the system can generate the 
observed flux efficiently. Our discussion is general and it does not depend on the specific energy 
generation or particle acceleration mechanism. 

We characterize the possible parameter space in terms of two parameters: the ratio of magnetic 
to electron's energy, e and the bulk Lorentz factor, V. A third auxiliary parameter is k, the ratio 
between the shell crossing time, and the angular timescale. Efficiency considerations alone limit 



the ratio of e: 10""^ < e < 10 (see ^3.5). The upper limit arises from the requirement that a 
large fraction of the total kinetic energy has to be stored at some stage in the relativistic electrons 
(eg > 0.1) and thus be available for radiation. The lower limit arises from the requirement that 
the up-scattered inverse Compton radiation does not dominate over the energy of the observed 
sub-MeV pulse. The lower limit becomes more stringent for shorter pulses and softer pulses (for 
example, tp ^ b x lO^^sec, Vp ^ lOOKeV lead to 8 x 10"^ < e < 10). Another implication of 
efficiency requirements is that soft - high luminosity bursts require huge intrinsic powers in the 
synchrotron model. This may explain why such bursts are not observed. 
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Efficiency also constrains T. The upper limit is determined by the condition that the emission 
process is efficient, namely that Uc < I'm- This is more constraining for small values of e and 
therefore it suggests rather strong magnetization. The lower limit on T arises from limits on the 
SSC component and from opacity considerations. These lead to 150 < F < 3000. Both the 
limits increase slightly with k, namely when the emitting radius is larger compared with the radius 
implied by variability. It is interesting to observe (see e.g. Figs. |3|4[ ) that within the allowed 
solution F is not necessarily correlated with e). 

Current Observations of high energy (GeV) photons lead to pair opacity constraints and to a 
larger lower limit on F. However, this limit depends on the fact that the high frequency photons 
originate at the same site as the sub-MeV photons (or at smaller radii) and it maybe relaxed if this 
is not the case. The observations of high energy photons actually only produces lower limits on 
the true frequency at which the system becomes optically thick for pair creation {vthick)- We find 
that if, in reality, this frequency is as high as SOOGeV, then low a bursts push the model further 
towards equipartition (e > 5 x 10^"^ for A; = 1 or e > 5 x 10~^ for k = 10). 

We explore the implications on the synchrotron model of future high energy observations by 
the CTA, in the tens of GeV to tens of TeV energy range. We find three general possibilities. Most 
trivially, a null detection by CTA does not yield significant constraints on the model. Second, a 
weak detection would imply that Uthick is at least of the order of tens of GeV, and that there is a 
weak SSC signal up to this frequency. This increases the lower limits on F and e. Depending on 
the actual flux measured, these limits can be as strong as e > 0.3 and F > 500. Third, a strong 
signal in the CTA band (above the extrapolation of the LAT signal) would be suggestive of an SSC 
peak. This could allow us to determine at least two of the three free parameters of our model (F, 
e, k) and effectively determine the full solution or possibly rule out the synchrotron model in case 
either of the first two parameters turns out to be beyond the limits discussed in this paper. 

Even though we don't have a full solution we can obtain a reasonable estimate of the expected 
range of synchrotron parameters needed to produce the observed prompt emission. The emission 
radius, R, is relatively large: lO^^cm < R < 10^''cm and it is independent of e. The magnetic 
field, B, spans almost four orders of magnitude: 3 Gauss < B < 10^ Gauss. The lines of constant 
B are almost parallel to the Um = line. Larger values of B lead to faster cooling so that the 
value of B is almost proportional to Xms, the ratio between the number of instantaneously emitting 
relativistic electrons and the overall number of relativistic electrons. The allowed range for 7^ is: 
3 X 10"^ < 7m < 10^. These large values of 7^ (compared with mp/m^) imply that, in many cases, 
as we discuss below, only a small fraction of the electrons (10^° < Ne < 10^^) are accelerated 
to relativistic velocities and participate in the synchrotron emission process. These large values of 
7m also suggest that any SSC peak would be at energies of at least a few hundred GeV. This is 
above the LAT band, naturally causing the SSC signal to be hard to observe with current detectors. 
Interestingly these parameter ranges, are very weakly dependent on k and they remain within this 



31 



range as long as we keep the observed quantities at their canonical values. 



Slightly deviating from the generic approach we present here, we examined the implications 
of the model if most of the energy is stored at some stage in the form of kinetic bulk motion. This 
would be naturally the case in the internal shocks model but also if the jet is initially Poynting flux 
dominated but it converts it energy to kinetic energy before the energy is radiated away. In this case 
the synchrotron solution requires that only a small fraction of the electrons in the flow are heated 
by the shocks and dissipate their energy during the prompt emission ( [Daigne fc Mochkovitch 1998 



Bosnjak et al. 2009[ ). However, we note that the distribution of electrons needed for this type of 



solution is different than what was achieved by recent PIC simulations, as it requires a gap between 
the energy of the non-accelerated electrons, and the minimal energy of the power-law component 



(see 10). We find that the fraction of accelerated electrons, ^ must satisfy ^ < 4 x 10 "^k^/"^ for a 
proton dominated flow (see Fig. ^ and ^ < 2 x IQ^^k^/'^ (Fig. 1 1 1 for a pair dominated flow. 

The main caveat of this model, that we did not address here, remains the "line of death" 
problem, regarding the low energy spectral slope (see also Daigne et al.||2011] ). A partial solution 



to this problem can be achieved if the system emits in a marginally fast cooling regime, where 
z/„ ^ u^. The high V values required by opacity considerations, push the model towards this 
regime. However, even the marginally fast solution, only steepens a up to —2/3 and it does not 
solve the entire problem. A different option is to have a steepening of the slope due to high 



self absorption frequency (Granot et al. 2000). This is not easily achieved, as the self absorption 



frequency is typically four orders of magnitude below the peak. 

Overall we conclude that there is reasonable range in the parameter (F and e) phase space 
for which a synchrotron solution that fits the peak frequency and the peak flux is possible. The 
main characteristics of this solution are large emission radii, relatively high magnetization and 
large values of electrons' Lorentz factors. These Lorentz factors naturally explain the lack of a 
SSC signal within the LAT band. In addition, the solution naturally explains the lack of soft - high 
luminosity bursts. Possible future detections of a SSC signal at higher frequencies, will allow to 
either obtain the full solution of the model, or rule it out altogether. 

We thank Jonathan Granot, Rodolfo Barniol Duran and Indrek Vurm for many helpful discus- 
sions. The research was supported by an ERC advanced research grant (GRB). 
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